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We present an approach based on a dimer expansion which describes low-energy singlet excitations 
(singlons) in spin-1 Heisenberg antiferromagnet on simple square lattice. An operator (“effective 
Hamiltonian”) is constructed whose eigenvalues give the singlon spectrum. The “effective Hamil¬ 
tonian” looks like a Hamiltonian of a spin-- magnet in strong external magnetic field and it has 
a gapped spectrum. It is found that singlet states lie above triplet ones (magnons) in the whole 
Brillouin zone except in the vicinity of the point (7r,0), where their energies are slightly smaller. 

Based on this finding, we suggest that a magnon decay is possible near (tv, 0 ) into another magnon 
and a singlon which may contribute to the dip of the magnon spectrum near (tv, 0) and reduce 
the magnon lifetime. It is pointed out that the singlon-magnon continuum may contribute to the 
continuum of excitations observed recently near ( 7 r , 0 ). 

PACS numbers: 75.10.Jm, 75.10.-b 


I. INTRODUCTION 

Spin-1 Heisenberg antiferromagnet (HAF) on square lattice is one of the most extensively discussed models of 
quantum magnetism. Interest in this model is particularly stimulated by its relevance to the parent compounds of 
the high temperature cuprate superconductors* 7 Despite its simplicity and much theoretical and experimental efforts 
in studying this model and related compounds, it continues to give surprises. 

We discuss the simplest variant of this model which Hamiltonian has the form 

H = (1) 

(hi) 

where ( i, j) denote nearest neighbor spins and the exchange coupling constant is taken to be equal to unity. As HAF 
CD is defined on bipartite lattice, its ground state is a singlet according to the Marshall’s theorem*? - The singlet 
nature of the ground state does not conflict with the long-range Neel magnetic order obtained both theoretically and 
experimentally*! Although no direction is selected in the singlet ground state, the spin-spin correlation function is 
finite at large distances implying a net staggered magnetization^ The spontaneous symmetry breaking takes place 
in the strict thermodynamic limit accompanied by formation of Goldstone elementary excitations (magnons) carrying 
spin 1. 1 ' —— Semiclassical approaches based on spin-wave theory (l/S'-expansion) give surprisingly accurate analytical 
description of the majority of ground-state and low-temperature properties of this quantum models In particular, the 
low-energy part of the magnon spectrum £k is successfully described within the first order in 1 /S whereas contributions 
from higher-order terms are very small*!>* - — 

The focus of research is currently on short-wavelength magnons. Recent theoretical and experimental results show 
that £k is not reproduced quantitatively by analytical methods around k = (7r, 0). A roton-like dip is found in the 
magnon spectrum at (tv, 0) by Quantum Monte-Carlo,' 1 ^^ series expansion around the Ising limit*!! and continuous 
similarity transformation?! techniques (see Fig. [[]). Remarkable agreement between these very different numerical 
methods signifies that the anomaly in the spectrum around (tv, 0) is an intrinsic property of the model. Magnon 
spectrum extracted from the neutron scattering data obtained in the metal-organic compound Cu(DCOO )2 ■ 4 D 2 O 
(CFTD) which is a perfect realization of model (JT]) shows the roton-like dip at (7r, 0). This dip is quantitatively 
described by the numerical results mentioned above. 17-19 Apart from the magnon peak at energy transfer ui « 
£(^ 0 ) r: 2.19, a continuum of excitations is observed extending from the energy slightly below £(^, 0 ) (from lu r; 1.9) up 
towR 3.8*?! It is concluded in Ref*!! that some correlations lead to this continuum which are isotropic in spin space. 
Continuums of magnetic excitations around (7r, 0) have been observed also experimentally in some cuprates*?-? - ™ 
Fractional elementary excitations producing continuums in spin-spin correlators are frequent in quantum physics of 
one dimension*? 4 To account for the continuum of excitations, fractional spin-^ excitations (spinons) are discussed 
now in HAF as wel h 18 i 25 i 26 It is proposed that a magnon is a confined state of two spinons in the whole Brillouin 
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FIG. 1: (Color online.) Spectrum of spin-| HAF along high-symmetry paths of the Brillouin zone shown in the inset. Spectra 
of triplet excitations (magnons) are presented obtained in the third order in 1/ S^— by series expansion around the Ising limit pH 
and by Quantum Monte-Carlo (QMC) computation^ (available only for k = (n, 0) and (n/2,n/2)). Series expansion and 
QMC results are consistent with those found recently by continuous similarity transformation technique^ and they describe 
quantitatively the magnon spectrum observed experimentall y 17 ' 18 in Cu(DCOO )2 • 4 D 2 O. The spectrum is also presented of 
low-energy singlet excitations (singlons) obtained in the present paper. 


zone except in the neighborhood of the point (7r,0), where spinons are deconfined and they form the continuum 
of excitations. However, it is proposed in RefJ£ that all features discussed can be explained in terms of magnons 
interaction not invoking spinons (while this paper does not disprove spinons existence). 

Motivated by the observation of the isotropic excitation continuum which extends to quite large energy, we address 
the problem of low-energy singlet excitations in model ©■ To the best of our knowledge, only triplet excitations 
are discussed in the literature (spin-i spinons arise as parts of magnons) whereas a little is known about low-energy 
singlet excitations. We perform a sort of dimer expansion (“plaquette expansion”) to derive an operator (an “effective 
Hamiltonian”) whose eigenvalues give the spectrum of singlet excitations (singlons). Our starting point is a set of 
isolated plaquettes in which exchange coupling constants are equal to unity between all four spins (see Fig. [2j. To 
come from decoupled plaquettes to the square lattice, we introduce an operator controlled by a single parameter 
A that switches on interactions between spins from different plaquettes and weakens interactions between diagonal 
spins in each plaquette. Decoupled plaquettes and model 0 correspond to A = 0 and A = 1, respectively. Then, 
we perform perturbation calculations up to 7-th order in A and extrapolate results to A = 1 by means of standard 
methods developed for critical phenomena. We present arguments that it is the low-energy singlet sector of HAF CD 
that we consider in this way. 

It is clear that a quantum phase transition (QPT) occurs on the way from decoupled plaquettes to the square 
lattice because ground state is disordered and magnetically ordered at the beginning and at the end of this way, 
respectively. The ground state and low-energy excitations change drastically upon QPT that is accompanied by non- 
analytic behavior of physical quantities as functions of driving parameter near QPT.- 27 As the long-range magnetic 
order cannot arise in the first few orders of perturbation theory^ 8 dimer expansions are inappropriate for discussion 
of ground-state properties in the ordered phase. It is the reason why authors of previous dimer expansion s 28 ' 29 do not 
consider the ordered phase focusing on properties of QPTs from disordered phases to the ordered one. However singlet 
excitations lie well above long-wavelength triplet ones (see below) which determine properties of the QPT. As high- 
energy excitations normally do not change drastically upon phase transitions, this gives promise that dimer expansions 
can describe singlet excitations. That is why we develop below the special dimer expansion which is convenient for 
discussion of the low-energy singlet excitations. This expansion differs from those suggested in Refs^^S in which 
starting points are isolated couples of spins. 

The main result of our study (i.e., the singlon spectrum) is illustrated by Fig. [T] We find that the “effective 
Hamiltonian” looks like a Hamiltonian of a spin-i magnet in strong external magnetic field and that the singlon 
spectrum is gapped (the large gap value is in agreement with previous findings according to which the low-energy 
physics of HAF is governed by low-energy magnons). Singlet states lie above triplet ones (magnons) in the whole 
Brillouin zone except in the vicinity of the point (n, 0). Remarkably, the singlon spectrum crosses the magnon one 
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FIG. 2: Simple square lattice (a) and decoupled plaquettes having doubly degenerate singlet ground state (b). To come from 
decoupled plaquettes to the square lattice, we introduce operator © controlled by parameter A so that A = 1 and A = 0 
correspond to panels (a) and (b), respectively. 
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FIG. 3: Wave functions of an isolated plaquette from which its singlet ground state wave functions © are constructed. Bold 

lines denote singlet states of the corresponding two spins. 




at those points, where analytical results of the 1/5-expansion start to deviate from numerical data for the magnon 
spectrum. Based on this finding, we suggest that a magnon decay is possible near (n, 0) into another magnon and a 
singlon which may contribute to the dip of the magnon spectrum near (ir, 0) and reduce the magnon lifetime. Notice 
that such a microscopic mechanism fundamentally cannot arise in 1/5-expansion. We point out that one-singlon 
states are invisible for experimental methods measuring two-spin correlators. However triplet states containing one 
magnon and one singlon and forming a singlon-magnon continuum may contribute to the experimentally observed 
continuum of excitations near ( 7 r, 0). It should be stressed that our findings and conclusions cannot disprove neither of 
the physical pictures proposed before for the description of the high-energy peculiarities of model ©■ In the present 
stage, we can only point out the possibility of additional contributions from singlons to the discussed features. 

The rest of the present paper is organized as follows. We describe our approach in Section [TT] The “effective 
Hamiltonian” is constructed and analyzed in Sections m and IIV1 respectively. Results obtained are discussed in 
Section [V] Section IVT1 contains our conclusion. 


II. PLAQUETTE EXPANSION 

Wave functions of the doubly degenerate singlet ground state of isolated plaquette 

*+ = ^1 +^2 

’ ( 2 ) 

AR” = </>l — 4>2 

are constructed as linear combinations of nonorthogonal ones </> 1 = (| f)i| 4-)2 — | l)i| t) 2 )(| t)31 4-}4 — I -l)3| t} 4)/2 
and </>2 = (| t} 2 | l }3 — | D 2 I t) 3 )(I t) 4 1 |)i — | ^ 4 1 t) 1 )/2 which are depicted in Fig. [3] It is convenient to consider 
and as states of a pseudospin s = \ corresponding to s z = 1/2 and s z = —1/2, respectively. 

To come from decoupled plaquettes to the square lattice, one has to switch on the following spin interactions which 
magnitudes are controlled by the single parameter A: 

V = A (S^S« + S^si j) ) + (s^ )s i P) + S«S,«) - 

<*,P> 

where upper and lower indexes of S enumerate plaquettes and the spin number in the plaquette (according to Fig.©, 
respectively, and (i,j) and (i,p) denote nearest neighbor plaquettes in horizontal and vertical directions, correspond¬ 
ingly. Decoupled plaquettes and HAF on the square lattice m correspond to A = 0 and A = 1, respectively. 



(3) 
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The singlet ground state is 2 N times degenerate in the system containing N decoupled plaquettes. This degenerate 
energy level splits into a singlet band (at N —► oo) upon A increasing because V commutes with operator of the 
total spin. Our goal is to describe this band in terms of interaction between pseudospins. It is our hope that at 
least lower part of this band forms the low-energy singlet sector of HAF (JT]). It may happen, however, that some 
upper singlet levels cross all energy levels stemming from the lowest singlet one upon A increasing from 0 to 1. 
In this case, all the band of singlet excitations we consider below would describe some excited states of HAF. We 
suggest the following symmetry arguments against this scenario. It is important that ’F+ and tH - given by Eq. © 
belong to different irreducible representations of the point group C± v of the plaquette. Then, in any cluster having 
the symmetry C^ v and containing N plaquettes (see Fig. [IJb) for the smallest cluster of this kind), one can find 
functions belonging to any irreducible representation of the group C± v among 2 N ground state wave functions. If the 
perturbation is one-parametric, the level crossing is generally forbidden of two energy levels belonging to the same 
irreducible representation. 30 Then, the lower singlet levels stemming from the lower degenerate level of such clusters 
would repel all the upper singlet levels upon A increasing from 0 to 1. Bearing in mind that one can consider arbitrary 
large clusters of this kind and that boundary conditions in a large enough system have normally a very small impact 
on its bulk properties, it is reasonable to conclude that we deal with the singlet band at least lower part of which 
forms the lowest singlet sector of HAF ([!]) . 


III. “EFFECTIVE HAMILTONIAN” CONSTRUCTION 


Our main goal is to obtain the operator H (the “effective Hamiltonian”) acting in the pseudospin space which eigen¬ 
values give energies of levels in the lower singlet band (assuming that N —> oo). Consideration of the corresponding 
wave functions in the spin space is out of the scope of the present paper. We use the standard perturbation theory for 
systems with degenerate energy levels developed by Bloch and described, e.g., in the textbook^. The conventional 
series expansion technique has a similar background— 

To illustrate the idea, let us consider corrections of orders A 0 and A 1 . In the zeroth order, H = CN , where 
C = — A°3/2 is given by the ground state energy of isolated plaquette. The first order corrections to energies of lower 
singlet levels are determined by matrix elements (4'f4'^...4 , ^|V|'Ff4'^...d'^), where V is given by Eq. ([3]) and 
denote states © of i-tli plaquette. It is easy to check that the first two terms in Eq. ((3J transform a ground state 
wave function of a plaquette to the triplet sector. As a result nonzero corrections in the first order in A originate only 
from the third term in Eq. ([3]) which contains interactions between spins belonging to the same plaquette. In terms 
of pseudospins, these corrections can be described by an interaction of each pseudospin with an effective magnetic 
field h: H = CN + where C = -A°3/2 + i(('I'-4'-| V 3 |'p-'I'-) + ('F+T+|y 3 |'F+'F+}) = -A°3/2 + A/2, 

h = (4'~4 ,_ |V 3 |\k~4 ,_ ) — ('F+<F+|y 3 |vI/+vl/+) = 2A, and U 3 = — A(SiS 3 + S 2 S 4 ) (cf. the third term in Eq. @). 

Terms of the perturbation theory of the second order in A contribute to C, h, and they lead also to terms describing 
an interaction between nearest-neighbor pseudospins. As a result H acquires the form in the second order in A 


\W 


\ ' 


H = CN + h st + ^ 


+ s + 
n “* 3 


-\-J A. 


Si Si 


s i sZ j + Jtj~ s ts] + J i j Z s i Sj + J~+si + JA sfSj + J^ + St 

(i,j) 

+ Jii + Si S+ + J, 


T+ ~sfs~). 


(4) 


Apart from corrections to coefficients in Eq. 0. higher-order terms of the perturbation theory lead also to multi¬ 
pseudospin interactions and to two-pseudospin long-range interactions. 

We calculate corrections to all parameters in the “effective Hamiltonian” H (including the multi-pseudospin and the 
long-range terms) up to the 7-th order in A. Series are presented in Table Q] for all the nonzero coefficients excluding 
numerous coefficients for long-range and multi-pseudospin interactions. It should be noted that in accordance with 
the general property of the perturbation series,— operator H is non-Hermitian (see, e.g., series for J ++ and J ). 
We find below that the spectrum of H is real as it must be. Operator H is translationally invariant: it is defined on 
the square lattice which period twice as large as the period of the original lattice. 39 

To extrapolate the spectrum from A <C 1 to A = 1, we use Pade and Pade-Borel resummation techniques —— There 
is no point in applying these techniques individually to each coefficient in H because the number of terms is small in 
the series for some of them (see, e.g., Table HI . Then, we derive below analytical expressions for physical quantities 
and find series for them up to the 7-th order in A using series for the coefficients in H. 
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TABLE I: Nonzero coefficients of the “effective Hamiltonian” H in the first seven orders in A describing interaction between 
nearest and next-nearest neighbor pseudospins. Numerous coefficients for long-range and nrulti-pseudospin interactions have 
been also calculated but they are not presented here. Subscripts h, v, and d denote the shortest horizontal, vertical, and 
diagonal bonds, respectively. The constant C has also the zero-order correction equal to —A°3/2 which is given by the ground 
state energy of an isolated plaquette. 



l 

2 

3 

4 

5 

6 

7 

c 

0.5 

-0.416666 

-0.250000 

-0.266783 

-0.315197 

-0.486108 

-0.807964 

h 

2.0 

0.166667 

-0.375000 

-0.574073 

-1.031603 

-1.746714 

-3.619268 

JZZ 

J h 

0 

-0.083333 

-0.062500 

-0.218171 

-0.436133 

-0.937999 

-2.063095 

J h~ 

0 

-0.062500 

0.015625 

-0.031684 

0.031262 

-0.138602 

-0.195048 

A ++ 

0 

-0.062500 

-0.109375 

-0.219184 

-0.463558 

-1.208478 

-2.851198 

4 + - 

0 

-0.062500 

-0.046875 

-0.062934 

-0.058693 

-0.071379 

-0.063827 

Jh + 

0 

-0.062500 

-0.046875 

-0.062934 

-0.058693 

-0.071379 

-0.063827 

T~ z 

J h 

0 

-0.072169 

-0.018042 

-0.002506 

-0.027965 

-0.051093 

-0.238453 

J+ z 

J h 

0 

-0.072169 

-0.090211 

-0.110759 

-0.143319 

-0.206620 

-0.342623 

r~ 

j h 

0 

-0.072169 

-0.018042 

-0.002506 

-0.027965 

-0.051093 

-0.238453 

j z + 

J h 

0 

-0.072169 

-0.090211 

-0.110759 

-0.143319 

-0.206620 

-0.342623 

JZZ 

J d 

0 

0 

0 

0.012732 

0.054546 

0.010743 

-0.177531 

Jr 

0 

0 

0 

-0.009549 

0.002464 

-0.026063 

0.023463 

Jl + 

0 

0 

0 

-0.009549 

-0.040996 

-0.157527 

-0.502310 

ji + 

0 

0 

0 

-0.009549 

-0.019266 

-0.038716 

-0.058577 

Jl~ 

0 

0 

0 

-0.009549 

-0.019266 

-0.038716 

-0.058577 

JZZ 

J V 

0 

-0.083333 

-0.062500 

-0.218171 

-0.436133 

-0.937999 

-2.063095 

J v 

0 

-0.062500 

0.015625 

-0.031684 

0.031262 

-0.138602 

-0.195048 

Jt + 

0 

-0.062500 

-0.109375 

-0.219184 

-0.463558 

-1.208478 

-2.851198 

J H— 

O v 

0 

-0.062500 

-0.046875 

-0.062934 

-0.058693 

-0.071379 

-0.063827 

J ^ 

O v 

0 

-0.062500 

-0.046875 

-0.062934 

-0.058693 

-0.071379 

-0.063827 

]~ z 
a v 

0 

0.072169 

0.018042 

0.002506 

0.027965 

0.051093 

0.238453 

j+ z 

O V 

0 

0.072169 

0.090211 

0.110759 

0.143319 

0.206620 

0.342623 

r~ 

O v 

0 

0.072169 

0.018042 

0.002506 

0.027965 

0.051093 

0.238453 

jz + 
a v 

0 

0.072169 

0.090211 

0.110759 

0.143319 

0.206620 

0.342623 


IV. “EFFECTIVE HAMILTONIAN” ANALYSIS 

The resummation procedure gives for h quite a large value around 2 at A < 1. Then, it is reasonable to suppose 
that (s z ) ~ —1/2 in the “ground state” of H at A < 1 (this assumption is confirmed by the spectrum stability in 
further calculations). Then, it is convenient to use the conventional Holstein-Primakoff transformation to find the 
spectrum 


s i = y 2s — a\ai ai, sf = a\yj2s - a\a i: s* = -s + a\ai. (5) 

Substituting Eqs. © into H , one leads to the Bose-analog of the “effective Hamiltonian” H = E 0 H „, where 

Eq is the ground state energy not renormalized by pseudospin fluctuations and H n denote terms containing products 
of n operators a and a'. 

At least the first seven terms in the series for H\ are equal to zero. It is a consequence of the symmetry of series 
for coefficients contributing to Hi. This property can be illustrated by series for J ±z and J 2± presented in Table |U 
It is seen that series for these coefficients for vertical and horizontal bonds have opposite signs in all orders in A. One 
concludes also from this fact that quantum corrections to H i from terms H n with odd n vanish as well. 

The bilinear part of H has the form 

H 2 = Q k«-k + > 


( 6 ) 
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(c) 



FIG. 4: Simplest diagrams contributing to renormalization of the ground state energy (a) and to the singlon spectrum (b)-(c) 
due to fluctuations of pseudospins. 


where B m k ^ B* k so that H 2 is non-Hermitian. The Bose-analog of H can be analyzed in a standard way (see, 
e.g., Ref.—) by introducing Green’s functions G(k) = (ak,a k ) u , F(k) = (ak,a~k)w, G(k) = (al k ,a_k)w and Ft(fe) = 
(aL k , a k ) w , where k = (w, k). We have two sets of Dyson equations for them one of which has the form 


G(k) = G^(k) + G^(k)S(k)G(k) + G^(k)[B pk + n(k)]F*(k), 
Ft (fc) = G (0) (fc)E(fc)Ft (k) + G (0) ( k ) [B mk + nt (k)]G{k ), 


where G^°\k) = (ui — F k ) 1 is the bare Green’s function and E(fc), E(fc), II(fc), and n't (A:) are self-energy parts. One 
obtains solving Eqs. © and a similar set of equations for G(k) and F(k ) 


where 


G(k) 

F(k) 


ui + £k -(- E (k) ^ _ — ui + Fk + E(fc) 

V{k) ’ V(k) 

B p k + n(fc) _t/7\ F m k + nt(fc) 
v(k) ’ V{k) ’ 


V(k) = uj 2 ~ e 2 0k - fl(k), 

^ok — \j F k B mk B pk , 

D(fc) = F k (E + E) - F mk n - F pk nt - w(E - E) - IUlt + EE, 


( 8 ) 


(9) 

( 10 ) 

( 11 ) 


G[k) = G(—k), E (k) = E(— k), and eok is the spectrum in the linear spin-wave approximation. Quantity D(fc) given 
by Eq. m describes renormalization of the spectrum square. We find D(fc) within the first order in 1/s by calculating 
corresponding diagrams for self-energy parts shown in Figs, ©b) and[4](c). The diagram presented in Fig. |4](a) gives 
the first 1 /s correction to the ground state energy. All the diagrams give small contributions to the physical quantities 
due to large gap in the bare spectrum eok (see below). They are calculated as follows. We carry out integration over 
internal energy, expand the resultant expression into a series up to the 7-th order in A, and carry out the summation 
over internal momenta. The external energy is taken to be equal to eok in the diagram shown in Fig. [2c). Due to the 
smallness of 1/s-corrections, A is actually the only parameter in our theory controlling series for physical quantities. 


V. RESULTS AND DISCUSSION 

Ground state energy. The series for the ground state energy per spin E gs not renormalized by pseudospin fluctuations 
is obtained from H by replacing s\ and sf by —1/2 and 0, respectively. Then, the constant term and terms containing 
only ^-components of pseudospins contribute to E gs . Correction from the diagram shown in Fig. ©a) leads to only 
slight variation of coefficients in the series (that does not practically change the result of resummation) which has 
the form -0.375 - 0.125A - 0.135417A 2 - 0.0239258A 3 - 0.0262953A 4 - 0.00190434A 5 - 0.0144457A 6 + 0.0155128A 7 
(all coefficients in all series are calculated with machine precision). As coefficients in this series do not rise and have 
irregular signs, we analyze only Pade approximants and obtain E gs = —0.694(5). The estimated errors in all the 
extrapolated values are always somewhat subjective— We exclude first all the approximants with obvious “defects” 
containing spurious poles with low residue at small A. The remaining approximants are weighted to favor the higher- 
order ones. The error is estimated from the scatter among the weighted values. 

Previous numerical studies^ ' 16 i 36 i 37 give values for the ground state energy around —0.669 that differs slightly 
from our finding. We attribute this discrepancy to the above discussed inability of dimer expansions to describe the 
ground-state properties in ordered phases. As the long-range magnetic order cannot arise in the first few orders of the 
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q 

FIG. 5: Illustration of the possible microscopic mechanism contributing to the dip in the magnon spectrum near k = (7r,0): 
magnon decay into another magnon and a singlon. Solid and wavy lines stand for magnon and singlon Green’s functions, 
respectively. 


perturbation theory, one has to analyze the whole series in A to find, e.g., the ground state energy at A = 1. Other 
dimer expansions underestimate E gs as well in the ordered phase of HAF *22, 

It should be pointed out also that contributions of long-range and multi-pseudospin interactions are not negligible: 
the value of E gs found using “effective Hamiltonian” © with coefficients from Table IHis -0.685(5). 

Spectrum of singlet excitations. The spectrum is presented in Fig. [I] found from the bilinear part of the “effective 
Hamiltonian” © (i.e., using Eq. (1101) 1 by means of Pade approximants analysis as well as Pade-Borel resummation 
approach. In particular, series for the spectrum at k = (7r, 0) and (7t/2, 7t/2) have the form 

e (jr , 0 ) = 2A + 0.083334A 2 - 0.453125A 3 - 0.356698A 4 - 0.47287A 5 - 0.342342A 6 - 0.468077A 7 , (12) 

£( 7 r/ 2 , 7 r/ 2 ) = 2A + 0.583334A 2 - 0.078125A 3 + 0.254847A 4 + 0.0564162A 5 + 0.253933A 6 - 0.173274A 7 (13) 

which give after the resummation 1.76 ± 0.15 and 2.75 ± 0.15, respectively. Again, contributions of long-range and 
multi-pseudospin interactions to these results are not negligible being of the order of 10%. As the gap in the singlet 
spectrum is large, one expects a small spectrum renormalization by pseudospin fluctuations. For instance, diagrams 
shown in Figs.QJb) and©c) do not practically change the results. 

It is interesting to speculate about the role of singlet excitations observed above using the quasiparticle concept 
thus assuming that each excited singlet state corresponds to a quasiparticle (singlon) carrying spin 0. It should be 
noted that one-singlon states are invisible for experimental methods measuring two-spin correlators. For instance, 
they give zero contribution to the dynamic structure factor oc Im(<Sk, S^_ k ) u because any operator 5“ moves the 
singlet ground state wave function to the triplet sector. However triplet excited states containing one magnon and 
one singlon can contribute to two-spin correlators. Such states can form a continuum at any given momentum k 
determined by the sum e q + ek^ q which starts at the singlon energy tk and ends at £k + £o- Using our results for ek 
and thosei2 for the magnon spectrum £k, we obtain at k = (7r,0) that this continuum lie within the energy interval 
(1.76 ± 0.15, 3.95 ± 0.17) which is in agreement with that (~ 1.9, ~ 3.8) observed in CFTD experimentally^ A weak 
continuum is seen also in the experimental data at (7 t/2 ,7 t/2 ) lying in the range (~ 2.3,^ 3.8). Our results give the 
interval (2.75 ± 0.15,4.14 ± 0.17) in this case. Thus, the singlon-magnon continuum can contribute to that obtained 
experimentally. 

Curiously, the singlon spectrum shown in Fig. [I] lies below the magnon one at that region of the Brillouin zone, 
where analytical results of the 1 /iS-expansion deviate from numerical data for the magnon spectrum. Then, one can 
speculate about the microscopic mechanism contributing to the dip in the magnon spectrum near (7T,0). As the 
magnon spectrum enters into the singlon-magnon continuum at that region, a magnon decay may become possible 
there into another magnon and a singlon (see Fig. ©. In many cases, such decay processes lower the quasiparticle 
energy and reduce its lifetime. Notice also that such a microscopic mechanism fundamentally cannot arise in 1/S- 
expansion. Derivation of an effective Hamiltonian describing both singlet and triplet excitations of the initial model 
© would be an interesting subject of further discussion. Such effective Hamiltonian would provide a quantitative 
description of the processes shown in Fig. [5] 

To conclude, singlet excitations may contribute to the high-energy peculiarities of HAF © and further discussion 
is required to fully clarify their role in this model. 


VI. CONCLUSION 

We develop the approach based on the plaquette expansion to discuss the spectrum of lower singlet excitations 
in spin-i Heisenberg antiferromagnet © on simple square lattice. The operator is found in the first seven orders 
of the perturbation expansion whose eigenvalues give the spectrum of low-lying singlet excitations (singlons). The 
resummation procedure gives the singlon spectrum shown in Fig. [1] We obtain that lower singlet excitations lie below 
triplet ones near the point k = ( 7 r, 0). We suggest that the magnon decay is possible near (7r, 0) into another magnon 
and a singlon (see Fig. [5]) which may contribute to the dip of the magnon spectrum near (n, 0) and reduce the magnon 




lifetime. This microscopic mechanism fundamentally cannot arise in l/S'-expansion. It is pointed out that one-singlon 
states are invisible for experimental methods measuring two-spin correlators. However the singlon-magnon continuum 
which can arise in this model may contribute to the continuum of excitations observed recently in CFTD. 18 
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